Correlation functions in decorated lattice models 
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Occupation probabilities for primary-secondary-primary cell strings and correlation functions for 
primary sites of a decorated lattice model are expressed through the well-studied partition function 
and correlation functions of the Ising model. The results are analogous to those found in related 
lattice models of hydrophobic interactions and are interpreted in similar terms. 

PACS numbers; 05.70.Ce, 05.70.Fh, 61.20.Gy. 



I. INTRODUCTION 
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Decorated-lattice-gas models that are equivalent to an underlying Ising model have been important as models of two- 
component mixtures with closed-loop solubility curves . In summing over the states of the decoration sites in the 
partition function one obtains an Ising model in which the energy and field parameters are related to the parameters 
of the original mixture model by certain rules (transcription relations). The interactions in the Ising model or in the 
equivalent one-component lattice gas may then be understood to have arisen, or to have been altered, through the 
mediation of the particles occupying the decorated sites (cells). In particular, the correlation functions between the 
particles occupying the primary cells, which would be those of the underlying Ising model or one-component lattice 
gas, may be understood to have resulted from such mediation. This is then analogous to the circumstance in recent 
lattice models of hydrophobic interactions , in which the solvent-mediated potential of mean force between solute 
molecules is obtained in terms of the correlation functions of the pure solvent. In the present work the potentials of 
mean force between occupants of the primary cells in a class of decorated-lattice-gas models with closed-loop solubility 
curves are calculated and interpreted in terms similar to those in the lattice models of hydrophobic interaction 



The model and the calculational machinery arc outlined in the next section and the results are then illustrated by 
numerical examples in Section III. In Section IV the correlation lengths, which are the exponential decay lengths of the 
various correlation functions, are calculated from those of the underlying Ising model and are displayed numerically. 

II. DECORATED-LATTICE-GAS MODEL 

We start with a brief review of a version of Mermin's decorated lattice model , Fig. 1, that describes a liquid mixture 
that possesses upper and lower critical solution points and a closed-loop temperature-concentration coexistence curve 
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Fig. 1. Mermin's decorated lattice model 
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Each cell, primary or secondary, is occupied by a molecule, either of type 1 or 2. Each molecule of either type has 
uj possible orientations. The only interaction takes place between the molecules in adjacent primary-secondary cells. 
This interaction is defined to be unless adjacent primary-secondary cells are occupied by molecules of different types, 
and its value then depends only on the orientation of the molecule in the secondary cell. For each primary-secondary 
neighboring pair, if the occupant of the secondary cell points to the primary cell occupied by an unlike molecule, the 
energy of interaction is U2 < 0; if it points in any other of — 1 directions, the energy is Ui > 0. 

To calculate the partition function for the whole system we first write the partition functions Qijk for all primary- 
secondary-primary cell strings for fixed occupations i and k of primary cells, 

Qui = Q222 = ^ 

Q121 = Q212 = 2e — — + (cj - 2)e-H- (1) 
Q122 = Q211 = Q112 = Q221 = e kT -I- (w - l)e kT , 



We also introduce the activity ratio C — Zi/z2 



mi d 

C = e [ ]2, (2) 

TO2 



where rrii and fii are molecular masses and chemical potentials for each species and d is the dimensionality of space. 
S in Eq. (62) of Ref. |^], we express the partition function of the model for the fixed number N of primary cells as 

y = u;^^C'^^(gniC + 0i2i)'^" (3) 

(g212C + Q222)'^^^(Qll2C + ^122)"^^^ 

Here Ni is the number of molecules of component ?, Nij is the number of i,j pairs of neighboring primary cells, and 
the sum is over all possible occupations "{rii}" of the primary cells. 

Now we can evaluate the probability Pijk to find a primary-secondary-primary cell string being occupied by the 
molecules i, j, and k. It can be written as 

(QiiiC + Qi2i)^'' (Q212C + Q222)'^^^(Qii2C + ^122)"^^^ 

Here ^ is the total number of bonds between N primary cells of the lattice with coordination number q. If we 
multiply Qijk by an auxiliary factor ayfe, so that Qijk is replaced by QijkCKijk, the probability Pijk can be expressed 



^ _ 2 dliiY 
qN daijk 



(5) 



We also can calculate the number of occupants of secondary cells pointing towards the opposite species occupants 
of primary cells, Np. This quantity is a measure of orientational ordering of the system. Note that the orientations of 
molecules in primary cells do not enter the expression for the partition function and therefore are isotropic everywhere 
in the phase diagram. Since every secondary-primary cell neighboring pair where the cells are occupied by opposite 

-[/2 

species and secondary cell occupants point to the primary cell comes with a factor , there is no need to introduce 
any auxiliary multipliers: 

""^--kf^ 

The structure of this expression is similar to that of it gives the probability of finding a particular ordered 
primary-secondary-primary string multiplied by the average number of such strings. 

To get the corresponding intensive quantity Up, we need to divide Np by the number of primary-secondary cell 
neighboring pairs qN (twice the number of primary-primary neighboring pairs), rip = Np/qN . 

Another important class of correlation functions it is possible to calculate is a connected pair and higher-order 
correlation functions for the primary cells. To underline the connection to the Ising model, let us introduce the 
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primary cell occupation numbers, Uj = +1 for the jth cell being occupied by a molecule of species 1 and rij = — 1 if 
it is occupied by a molecule of species 2. Then A^i and N2 in Eq. ^ can be expressed as Ni = {N + J2j=i and 

N2 = {N — We also introduce an auxiliary local field 7j, in the presence of which the partition function 

takes the form: 

r = w^(;^/2 ^ e5Er=i"'('^'+^i"C) (7) 
{".} 

aNii AN22 aNi2 
^11 ^22 ^12 • 

Here we used the shorthand notations: An = QmC + Q121, M2 = Q212C + Q222, M2 = Q112C + Qi22- 

As in the Ising model, the two-body connected correlation function for the primary cells i and j, C^j- = (jiiUj) — 
{ni){'nj), can be expressed as 



(8) 

{7fe}=0 



In order to complete the calculation of correlation functions Pijk and , we express the partition function i^({7j}) 
through the known partition and correlation functions for the Ising model. Using lattice identities, qNi = 2Nii + N12 
and qN2 2N22 + iVi2, we reduce Eq. (0) to 

r = ^^C^/2(AnA22)^ ^e^^"."^(-^'+^i"C+fi"4^) (9) 

{".} 

On the other hand, for the Ising model with the Hamiltonian 

H = -J^S^SJ-^S^(i^ + 7^), (10) 

the partition function Z and the connected two-point correlation function Cl^^""^ can be expressed as 

^e^S"i"-(-f^+'^')+2-^^+-, (11) 



_ qN_ 

e 



{S^} 



Here iV_| is the number of bonds between up and down spins, which is equivalent in the mixture model to the 

number of 1-2 neighboring pairs of primary cells, A^-i = Ni2. 

We observe that the connected correlation function for the decorated lattice model and Ising model are identical if 
the coupling constant J and external field H in the Ising model are defined as 

c,, = c^r' (13) 

*^<=>-5W^^], ^<=>ilnC+|ln^. 

We next proceed to the calculation of the /c-point potential of mean force between the occupants of the primary 
cells, Wk{{i ■ ■ . j}) — —kT ln[ ^'^^^'^'f^^ ]. Since the effective interaction between occupants of primary cells depends on 
the species and state of the occupants of the intervening secondary cells, the corresponding primary-primary cell 
correlation functions may be viewed as having been mediated by those secondary-cell occupants. This is analogous 
to the solvent-mediated correlations between solute molecules in related hydrophobic-interaction models |4|-||] , but it 
is to be understood that in the present model the effective correlations are not those between molecules of the same 
species since molecules of both species may be the occupants of either kind of site. Indeed, the correlations calculated 
here are those between the occupants (whatever their species) of the same kind of site (the primary cells). 

In the definition of W , pi is the non-negative occupation number (density) of the ith cell. Since the actual amplitude 
of p in the definition of W cancels, it is convenient to let it vary between and 1, i.e. to define it as pi = "'^^^ . 



(12) 
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It allows a fc-point potential of mean force to be expressed through the k- and lower-order connected correlation 
functions in the Ising model. For example, for k = 2, 



W2{{i,j})^~kT\n[^^£^ + l]. (14) 

Here (n) — ^g'^ is the average occupation number or magnetization for the corresponding I sing model. The corre- 
spondence between the parameters of the decorated lattice model and Ising model is given by (|l^). Similarly, one can 
calculate the higher-than-second order correlation functions and potentials of mean force. 

Finally we show how to calculate the primary-secondary cell and secondary-secondary cell correlation functions and 
potentials of mean force. Taking into account Eqs. (|l]), we calculate a probability S^njj to find a particle of species 
1 in the secondary cell situated between primary cells occupied by molecules of species i and j: 

S{iii} = -{Uj^+u^t ^2ur^ 

u:+2e kT +(uj~2)e kT 
S{112} = S{i22} = 5, (15) 

rt _ 2e kT +(u-2)e kT 

^{212} — -(C/1+U2) -^2UT - 

uj+2e kT 2)e kT 

Taking into account the definition of p for primary cells given above, we express an average density of molecules of 
the species 1 in secondary cells, p\ as: 

P'^ = {P'iPj)S{lll} + 2((1 - Pt)Pj}S[211} + ((1 - - Pj))S{212}- (16) 

Here pi and pj are the occupation numbers of the neighboring primary cells i and j. Similarly, the priniary-secondary- 
cell correlation function {pkp\) is expressed through the following average: 

{Pkp\) = {pkPiPj)S{iii} + {pk{l - Pi)Pj)S{2ii} + (17) 

(PfePi(l - Pj))S{ii2} + (Pfe(l - - Pj))S{2l2}- 

Here i and j are the primary cells between which the secondary cell / is situated. It follows that the primary- 
secondary cell correlation function of the decorated lattice model can be expressed through the energy-magnetization 
(3-point) correlation function of the underlying Ising model and various 2-point correlation functions. Finally, for the 
primary-secondary cell potential of mean force W2{{k^ Z}), one obtains: 

W2{{Kl}) = -kT\n^^^. (18) 

. A similar approach shows that the secondary-secondary cell correlation function and potential of mean force for 
the decorated lattice model can be expressed through energy-energy (4-point), energy-magnetization (3-point) and 
2-point correlation functions. 

Because of the generic nature of Mermin's decorated lattice model, these results are applicable for lattices in 
arbitrary dimension d with various coordination numbers q. 



III. EXAMPLES 



Using the formal mapping of the correlation functions of the decorated lattice model onto the correlation functions 
of an Ising spin system, we perform computations of the 2- and 3- body potential of mean force in two dimensions. 
The following values were chosen for the model parameters: Ui = —U2 (with C/i, unspecified, then a scale factor for 
the temperature), lo = 100, C — 1- The phase diagram for these values of the parameters is shown in Fig. 2. 
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Fig. 2. Closed loop temperature-concentration solubility curve. The dashed line marks the concentration (p) sa 0.08 

Despite the fact that the correlation functions for the Ising model are available analytically in the form of series 
expansions (see, for example, [Q), we found it more convenient to generate them each time directly using Monte Carlo 
simulations. We used Metropolis and Wolff algorithms |^ with the lattice size varying from 100 x 100 to 200 x 200. 
All the results presented below were averaged over the whole system and over 200 configurations; configuration were 
considered different if there were separated by about 10 flips for each spin in the Wolff algorithm. The correlation 
functions and potentials of mean force were measured along the lattice ("normal") unit vectors ((1, 0) and (0, 1)) and 
along the main diagonal (1, 1). We found that in general the results for "normal" and "diagonal" measurements are 
very similar. An example is in Figure 3, where the two-body potential of mean force is shown. The dashed curve is 
a result of the superposition of normal and diagonal plots, with normal and diagonal points occurring at integer and 
\/2x integer lattice separation. 

The potentials of mean force W2 and W3 were calculated for the temperatures corresponding to the widest point of 
the phase diagram ((p) « 0.05, kT/Ui « 0.44), and also for the two temperatures corresponding to the intersections of 
the dashed line and the coexistence curve in Fig. 2, at both of which there is the same minority species concentration 
(p) w 0.08 {kT/Ui « 0.384 and kT/Ui w 0.513). The effective Ising coupling constant is J/kT = 0.47 for the first 
point and J/kT = 0.45 for the next two. 
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Fig. 3. Plots of the two-body potential of mean force at kT/Ui = 0.44 (solid line) and kT/Ui = 0.38,0.51 (dashed line). 

The potentials of mean force for primary-secondary cell occupants is shown in Fig. 4. The original primary cell and 
two primary cells that surround the secondary cell lie on the same lattice vector; the distance between primary and 
secondary cell is assumed to be equal to half-sum of the distances between the original and two neighboring primary 
cells. 
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Fig. 4. Plots of the two-body potential of mean force at kT/Ui = 0.44 for primary-secondary cells (solid line) and primary- 
primary cells (dashed line). 

For the three-body potential, we chose the two simplest configurations: In the first one, particles 1, 2, and 3 lie on 
a straight line with the distance r between particles 1 and 2 being equal to the distance between particles 2 and 3. In 
the second configuration particles 1, 2, and 3 form a right triangle in which the legs 1,2 and 1,3 are of equal length r. 
We are mainly interested in the question of how well the effective 3-body interaction is approximated by the sum of 
the three pair interactions. 

In Fig. 5 we present a plot of the difference between the sum of the three pair potentials and the true 3-body 
potential for a linear configuration, AW^3(r) = 2W2{r) + W2{2r) — W3{r,r,2r). Since all potentials W2 and W3 are 
negative for any r, one can observe that the sum of the three pair interactions clearly overestimates (the absolute 
magnitude of) the true 3-body potential. However, unlike in some one-dimensional models of hydrophobic interactions 
1^, this overestimate is less (in absolute value) than the interaction between the furthermost particles, 1 and 3. To 
show that, a plot of the corresponding two-body potential W2(2r) is presented in the same Fig. 5. 
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Fig. 5. Difference between the sum of the three pair potentials and the true 3-body interaction, AW3(r) = 2^2(7') -I- W2{2r) — 
W3{r, r, 2r) (solid line), and the pair potential W2(2r) (dashed line), for the linear configuration of points and the temperatures 
kT/Ui = 0.384 or 0.513 

Similar results for the triangular configuration of points are presented in Fig. 6. Here again, as in the linear case, 
the sum of all pairwise interactions overestimates the true 3-body potential, but, as in the linear case, by less than 
the magnitude of the pairwise potential between two most remote points, W2{'\/2r). 
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Fig. 6. Difference between the sum of the three pair potentials and the true 3-body interaction, AW3(r) = 2W2{r) + 
W2{V2r) — W3{r,r, y'l-) (solid line), and the pair potential W2{V2r) (dashed line), for the triangular configuration of points 
and the temperatures kT/Ui — 0.384 or 0.513 



IV. CORRELATION LENGTH 



In this section the correlation lengths for the decorated lattice model are calculated from those of the underlying 
Ising model. Before we proceed to the calculation, we note that the correlation lengths for any pair of species (i.e., 
1-1, 1-2, 2-2) are all identical in this model. The exponential range of correlation for a pair of molecules of species a 
and P occupying the primary cells is defined by 



1/W = - 1^1 

\R.\^oo 



(19) 



where hap(R) is the pair correlation function for species a and f3. However, all the pair correlation functions /iii,ft.i2, 
and /i22 are expressed as 



(20) 



where kap = (2(5q/3 — l)/4papp with pa = Na/N. Note that C{R) is the two-body connected correlation function 

Cij defined in Eq. (^ , but now expressed as a function oi R = Ri — Rj instead of as a function of the primary cell 

addresses Since the factor kaf3 is independent of R, all the correlation lengths for any pair of species in this model 
are identical to the length defined by 



1/5 = - lim \R\-'^\n\C{R)\. 

\R\~*oo 

On the other hand, the correlation length for the Ising model is defined by 

l/<e^'''"f = - lim |i?rMn|C^^™f(i?)|. 

\R\^oc 



(21) 



(22) 



Therefore the correlation length for the decorated lattice model is calculated from ^^^^f through the transcription 
given by Eq. (|l§. 

For the Ising model an accurate expression for the correlation length at zero field below Tc is available as a function 
of exp(-2J//cT) 0. Thus we can calculate ^ for the decorated lattice model along the coexistence curve from the 
upper to the lower critical solution point. For the simple cubic lattice, the following expression is used: 



i = 4cosh"^ I 1 



e / 



2A', 



(23) 
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1 ^1 



where for the simple cubic lattice in the direction ep = (1, 0, 0), 
and for the same lattice in the direction ei = (1, 1, 0)/^/2, 



201a; 



14 



and f = 1 



T_ 
16' 



32 



and f 



(24) 



(25) 



(26) 



The variable x in terms of which A2 is expanded is to be understood as ^12/^^11^22 for the decorated lattice model, 
instead of exp{—2J/kT) for the Ising model. 

Here we show the numerical results for the correlation length for the decorated lattice model of the simple cubic 
lattice. The results are obtained for various values for the parameter lu at \U2\/Ui = 1. Shown in Fig. 7(a) is the 
behavior of the correlation length ^ in the direction eb along the coexistence curves as a function of temperature 
between the upper and lower critical solution points. For any given lu, the correlation length takes a minimum value 
at a temperature lower than the midpoint between the upper and lower critical temperatures. With increasing uj, 
the minimum of the curve becomes deeper and shifts toward the lower critical temperature. This means that a 
larger number of possible orientations uj gives rise to a correlation length that increases more rapidly as the system 
approaches the lower critical solution temperature. Similar results are obtained for the direction ei as shown in 
Fig. 7(b). 
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Fig. 7. Behavior of the correlation length along the coexistence curve for the decorated lattice-gas model for various values of 
UJ at \U2\/Ui = 1. The lattice type is the simple cubic and the lattice direction in which ^ is calculated is (a) eb = (1,0,0) and 
(b) ei = (1, l,0)/\/2. Temperature is in units of Ui/k. 
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V. SUMMARY 



We have calculated and illustrated the correlations between primary cells for a decorated-lattice-gas model that 

exhibits a closed-loop phase diagram. The thrcc-body potential of mean force was also calculated and compared with 
the sum of the three two-body interactions. The calculations are related to those in earlier models of hydrophobic 
interactions and are interpreted in similar terms. 
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